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INTRODUCTION 


The  method  used  here  is  that  of  obtaining  the  rectanf  coordinates  of 
three  trackers  with  respect  to  a  coordinate  system  whose  x  y  plane  is  tangent 
to  ar.  ellipsoidal  earth  at  the  location  of  one  of  the  trackers.  Attention  is  called 
to  the  fact  that  the  calculations  for  this  are  executed  prior  to  any  real  time 
computations.  After  this  is  accomplished  it  is  analytically  simple  to  solve  three 
distance  equations  in  three  unknowns  to  obtain  the  rectangular  coordinates  of  the 
missile.  Complete  procedures  are  given  in  the  text. 
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Given  the  geodetic  latitudes  and  the  longitudes  of  the  trackers  (L  X^), 

T„  (L  AJ,  and  T„  (L  ,  AJ. 
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Figure  1 


Referenced  to  an  orthogonal  right-handed  Cartesian  coordinate  system  with 
its  origin  at  the  center  of  the  earth,  its  positive  Y  axis  passing  through  the 
meridian  whose  longitude  is  X^  ISO*,  and  its  positive  Z  axis  passing  through  the 
north  pole,  the  trackers  will  have  the  coordinates: 


cos  sin  ^  -  X^) 
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where 


j  =  1,  2,  and  3, 
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V  =  a  (1  -  e  sin  L^) 
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a  =  semi  major  axis  of  the  meridional  ellipse. 
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b  =  semi  minor  axis  of  the  meridional  ellipse, 
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Vj  =  distance  from  T\  along  the  radius  of  curvature  extended  to  its 

point  of  intersection  with  the  minor  axis  of  the  meridional  ellipse. 
The  chordal  distances  T  T  ,  T  T  ,  and  T  T  can  be  obtained  using  Euclidian 
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distance  formulas. 

The  distance 


|0Q1j=Vi  sinLj-Zj. 


Translating  from  0  to  yields  tracker  coordinates 


(X.  ,  Y  ,  Z.  )  , 
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where 


X.  =  X. ,  Y.  =  Y., 
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Z.  =  (V,  sin  L  -  V.  sin  L  -  Z  *  Z.)  +  Z 
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Next  rotate  through  an  angle  of  90s  -  about  the  X^  axis  and  then 

translate  from  and  T  .  The  coordinates  with  respect  to  a  coordinate  system 
with  its  x  y  plane  tangent  the  et'rth  at  T  ,  its  positive  z  axis  vertically  upward, 
its  positive  y  axis  north,  and  its  positive  x  axis  east  will  then  be  given  by: 
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y  =  *  Y2  sin  Lj  *  Z  cos  L  , 
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z  =  -  Y„  cos  L,  +  Z„  sin  L, 
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xo  =  X  , 
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y  =*Y3  sin  L  +  Z  sin 
Q1 
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The  preceding  work  is  all  precomputational;  i.  e. ,  none  of  it  is  done  in  real 


time. 


To  obtain  the  coordinates  of  the  missile  M  write  the  distance  equations: 
2  2  2  2 


x  +  y  +  z  =  r  , 


2  2  2  2  2  2  2 
x  +  V  +  z  -  2x2  X  -  2y2  y  -  2z2  z  +  x2  +  y2  +  z2  =  r2  ,  (6) 
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x  +  y  +  z  -  2x3  x  -  2yg  y  -  2z3  z  +  x3  +  y3  +  z3  =  r2  .  (7) 
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Subtracting  (6)  from  (5)  yields 

2  2  2  2  2 
2x„  x  +  2y„  y  +  2zn  z  =  r,  -  r  +  x.  +  >*  +  z" 
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and  subtracting  (7)  from  (5)  yields 

2XS  x  +  2y3  y  +  2z3  2  =  r  I2  -  r,,2  -  x*  *  y*  *  z  .*  . 
Solving  (8)  and  (9)  for  x  and  y  gives 
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32  23  32  2  3 


**Vs'V*) 


where 


,  2  2  2  2  2 
f  =  r  -r  +  x„  +  y„  +  z„ 
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f  =  r  -r  +  x„  *  y„  +  z_ 

3  13  3  J3  3 


Substituting  the  right  members  of  (10)  and  (11)  for  x  and  y,  respectively,  in 
(5)  and  collecting  like  terms  yields: 

4  *(X2y3  "x3sr9)2  +  ^9Z3  "y3Z2)2  +  (X'<Z9  'X9Za)2*  *  1 
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+  4  ^(f2y3  ' W  ^2Z3  _y3Z2)  +  ‘f3X2  ”  W  (”3Z2  ~ VV  *  Z  {H} 

Solving  this  equation  for  z  and  rejecting  the  root  bearing  the  negative  sign 


yields: 
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Figure  2 
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where,  in  (14) 
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a  =  the  coefficient  of  z  , 
b  =  the  coefficient  of  z  , 


c  =  the  sum  of  all  ether  terms  in  the  left  member. 


Substituting  the  right  member  of  (15)  in  (10)  and  (11)  yields: 


a<Vs  'W + Vi  -  v2>  (‘b+  Jb  -4ac) 
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The  x,  y,  and  z  coordinates  of  the  missile  car.  now  be  computed  from  (16), 
(17)  and  (15),  respectively,  in  real  time.  It  should  be  noted  that  in  (15)  the  a  is 
entirely  precomputed,  while  the  f  and  the  f  in  the  b  and  c  are  computed  in  real 

M  O 

time  from  (12)  and  (13),  although  parts  of  these  expressions  are  precomputed. 
Specifically  any  expression  involving  only  one  or  more  of  the  parameters  x  ,  y  , 

M 
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To  obtain  the  Cartesian  coordinates  of  the  missile  with  respect  to  another 
point  on  the  earth’s  surface,  use  transformations  (42)  and  (44)  from  A.  Saastad 
and  C.  F.  Kettner,  Screening  and  Processing  of  Mark  I  Azusa  Data,  Convair 
Astronautics  Report  AZN-26-055,  17  November  1958. 


Reference:  Bomford,  B.  C.  ,  Geodesy,  Oxford  U.  Press.  1952. 


6 


